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SUMMARY 

A probabilistic model utilizing random material characteristics to predict damage evolu- 
tion in textile laminates is presented in this paper. Model is based on a division of each ply into 
two sublaminas consisting of cells. The probability of cell failure is calculated using stochastic 
function theory and maximal strain failure criterion. Three modes of failure, i.e. fiber breakage, 
matrix failure in transverse direction, as well as matrix or interface shear cracking are taken into 
account. Computed failure probabilities are utilized in reducing cell stiffness based on the mes- 
ovolume concept. A numerical algorithm is developed predicting the damage evolution and de- 
formation history of textile laminates. Effect of scatter of fiber orientation on cell properties is 
discussed. Weave influence on damage accumulation is illustrated with the help of an example 
of a Kevlar/epoxy laminate. 



INTRODUCTION 

Textile reinforced composites represent an important class of advanced composite materi- 
als. Although their structural performance under uniaxial loading yields to that of unidirection- 
al composites, better impact resistance and balanced in-plane properties of these materials as 
well as ease of application and the resulting low manufacturing cost make them superior for a 
variety of applications [1]. Ishikava and Chou [1-3] have developed several micromechanical 
models for textile laminates. Unidirectional mosaic and fiber undulation models are used for 
effective thermoelastic properties calculation of plain weave fabric composites. The bridging 
model is developed to simulate the load transfer between the interlaced regions in satin weave 
composites. Although the elements of stress-strain analysis can be incorporated in the models 
and effects like “knee behaviour” of textile composites can be described [2], these models do 
not allow to analyze the damage evolution in textile laminates. 

It is well-known that texdle composites under loading undergo multiple damage formation 
on several structural levels that corresponds to their complex structural hierarchy. The process 
of damage accumulation in these materials is of stochastic nature due to a number of random 
factors, such as inevitable scatters in stiffness and strength properties of components, local 
structural inhomogeneities, and random loading. A combined effect of the above mentioned 
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factors causes scatter in static and dynamic strength of finite size specimens. On the other hand, 
for large specimens, where all the combinations of random factors are present, their effect may 
result in gradual accumulation of damage and a definite failure strength. In general, failure 
analysis of textile composites should involve consideration of several consecutive geometrical 
models, corresponding to micro-, meso-, and macrolevels of composite structure [4], 

Numerous scientists have studied composite strength problems with the help of statistical 
modeling (see [5-7] and references therein). Ovchinskii [5] proposed and worked out the funda- 
mental principles, models, and algorithms for stochastic simulation of fracture processes in 
reinforced composites, starting with micro-level consideration. However, the simulation proce- 
dures with micro-level consideration require extensive computational resources for laminate 
analysis. 

A statistical method for reliability analysis of laminates was developed in [8- 10], In these pa- 
pers, theory of random scalar and vector field excursions was applied for computing the prob- 
ability of excursions of a random stress-strain field beyond the limiting surface. The method 
was applied in [1 1-13] for ply-by-ply stochastic failure analysis and reliability calculation of 
laminated composite plates and shells. In these papers the whole composite layer was consid- 
ered as a basic element and maximum probability criterion was used for determination of dam- 
age sequence. 

Another approach for damage analysis of laminates was developed in [14,15]. The model 
is based on a division of each layer into statistically large number of mesovolumes. Concentra- 
tion of broken mesovolumes in plies is calculated as a probability of ply random strains to ex- 
ceed the limiting failure strains. These concentrations are utilized in ply stiffness reduction.The 
model allows to predict the gradual damage accumulation and deformation history of compos- 
ite laminates. 

At present, paper probabilistic modeling of damage evolution in textile laminates is per- 
formed based on the approach developed in [14,15]. 

PROBLEM FORMULATION 

Consider a textile reinforced laminated composite sketched in Figure 1, a. Each ply of 
this laminate consists of two sublaminas, made of a set of cells of two types (Figure l,b). The 
orientation of fibers in cell with respect to laminate axes, referred further as cell orientation, is 
described by angle 0. The relative orientation of cells within the same lamina is not necessarily 
perpendicular, in general. The weave type in this model is described by fraction of cells of par- 
ticular orientation in sublamina v. The fraction of cells of the other orientation is, obviously, 
l - v. Material properties of cells of different orientations within the same lamina can also be 
different, in general. 

Now, the problem can be formulated as follows. Assume that initial elastic material prop- 

~ k ~ k - k ~ % 

ertiesof cells are random E ,, E2 0 , Gn 0 , v 12 . Index k is the cell identification number. Hereandfur- 

ther tilda sign upper bar -, and D denote random characteristic, mathematical mean, and 
dispersion of characteristic, respectively. Load increment applied to such a laminate results in 
random stress-strain field in each ply, sublamina, and, consequently, in each cell. Even at very 
low levels of applied load a nonzero probability of cell failure exists and damages start to ac- 
cumulate in composite. Accumulation of damages causes reduction in laminate stiffness and 
redistribution of stresses between cells and plies. 
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figure i. Schematic representation of textile laminate (a) and lamina subdivision (b) 







Assume that in-plane stresses applied to composite are monotonically increasing func- 
tions of a parameter t. a, (r) = [o,, (/) ,o 22 (r) ,x 12 (0 1 , with parametric derivatives a,(t ) . Ran- 
dom deformations of composite can then be calculated using integral equations 

t 

e,(0 = py(t)a ; (x)i/t a) 

0 

where S :J are current effective laminate compliances. They depend on cell current elastic 
properties, ply weave and composite lay-up. 

Sij(x) =■■ L[£j(T),£*(x),c! 2 (x),v* 2 (t),e*,v*] (2) 

The current elastic constants are functions of the initial elastic constants and current dam- 
age functions in cells. 

e!(x),£S(x),g! 2 (x),vI 2 (x) = 

(3) 

r- * ~ * ~k k "I 

M[Ei 0 , £ 2 0 , Gi 2 0 , v, 2 o , r i (x) J 

Damage functions may be calculated using cell random stress-strain field parameters and 
some appropriate failure criteria. 


r*(x) = £[e*(x),o*(x)] (4) 

Cell stress-strain field parameters are calculated, in turn, from known composite strains e, 
(eq. 1). 


Mx),o* (X) 


*[£,. (x) , E x (x) , E 2 (x) , G n (x) , v n (x) , 9*. v*] 


(S) 


L, M, R, K, P are stochastic functional operators to be specified. According to this ap- 
proach, current composite elastic properties, and therefore composite deformations and dam- 
age functions, are dependent on loading history. In order to integrate eq.(l), we have to 
calculate stochastic stress-strain field parameters in cells which depend on stochastic material 
properties and loading, as well as deformation history. The failure criterion (eq. 4) should be 
chosen to obtain probabilities of ply failure (damage functions). The stiffness reduction algo- 
rithm due to damage accumulation in plies (eq. 3) has to be specified. 


FAILURE PROBABILITY CALCULATION 

-k -k -k _* 

Assume, that cell stochastic elastic characteristics £,,£ 2 , Gi 2 ,v 12 are statistically indepen- 
dent and normally distributed random numbers. They can be described using mathematical 
mean values e\, e\, Gp, v* and dispersions D ., D ., D . , D , . Mathematical mean value and 

dispersion of a function of independent random variables y = f(x u x 2 ,-,x„) can be evaluated as 
follows: 
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( 6 ) 


y=f(xvx *«) 



Stochastic cell strain response to the incremental laminate loading is calculated using nor- 
mally distributed numbers properties. Mean values and dispersions of the following random 
variables are computed in sequence (see [14,15] for some of the formulas involved): 

(i) cell current technical elastic characteristics 

(ii) cell stiffness coefficients in material axes 

(iii) cell stiffness coefficients in laminate axes 

(iv) sublaminas stiffness coefficients 

(v) laminate stiffness coefficients 

(vi) laminate compliance coefficients 

(vii) laminate strain increments 

(viii) total laminate strains 

(ix) sublaminas stresses 

(x) cell strains in laminate axes 

(xi) cell strains in material axes 

Cell properties averaging procedure according to weave parameter v (iv), and classical 
lamination theory are used in the calculation. Procedure for current elastic properties estima- 
tion (i) is described below. 

Maximal strain criterion is applied to calculate the failure probability of individual cells. 
We will describe damage criterion for deterministic and stochastic failure strains. In case of de- 
terministic failure strains, each cell is assigned limiting failure strains in tension and compres- 
sion, e'*, e" *, respectively. These limiting values are the same in magnitude for shear strain. In 

case of stochastic failure strains these are normal distributions. The probability of i* to exceed 
deterministic border e’* is simply the area under differential distribution of e, on the right side 
of the upper bound e'*, and can be calculated from integral distribution F £ , of ply strain by using 
eq.(7). 


r*= 1 - F t (e'*) 


( 7 ) 


- k 


For stochastic failure strains e'- , probability i* to exceed e’ ( . can be calculated as probabil- 


- k 


. k 


ity of random difference e', - £, to be negative. 



Parameters of integral distribution F , » can be calculated as 


(8) 



D 

E 


,k 



+ D 


( 9 ) 
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Taking into account both upper and lower bounds e'* and e"‘, and assuming that strain field 
excursions beyond the two bounds are statistically independent events, the following expres- 
sions for rf can be derived 

(a) deterministic bounds e'*, e" * 


r- = 1 -F >&'*)+ FAST-) 


~ k - k 

(b) random bounds e' ; , e" , 


(i«) 


= F 


.» k (0) + F t ..,(0) 
£ - £ £*-£ ‘ 


(ID 


It is easy to show that second case is reduced to first one when D , t =£>„* = 0. Thus, dam- 

age state of each cell is described by three damage parameters r* = { A, / 12 } . Note, that ac- 
cording to the developed model damage state in cells depends on weave of textile reinforcement 
and for v * 0.5 , cells with the same orientation in top and bottom sublaminas, which have iden- 
tical initial elastic characteristics, may become fairly distinct due to damage. Hence, each ply 
of textile laminate in the proposed model is composed of cells of four different types. 


EFFECT OF DAMAGE ON CELL PROPERTIES: MESOVOLUME CONCEPT 

In order to use the information about the failure probability we need to specify how is it 
realized in the composite, i. e. when and how much damages of different types do emerge and 
how do they effect the current mechanical properties of composite. Assume that there are three 
types of damages related to three in-plane deformations in each cell e* = [e*,, e 22> y^] . Failure 

type due to e*, exceeding the critical bound can be associated with fiber breakage, and those 
due to e 22 and yf 2 - with matrix cracking in transverse to fiber direction and in shear, respective- 
ly. 


Assume that each cell consists of certain number of mesovolumes. The mesovolume has 
to be structurally homogeneous, i.e. to contain sufficient amount of reinforcing fibers to be 
considered as a continuum, and at the same time it has to be comparatively small in order to 
satisfy the condition of stochastic homogeneity of stress and strain fields inside its envelope. 
Initial random characteristics of cells are equal to respective mesovolumes characteristics. As- 
sume that each of the mesovolumes can be either perfect or broken in any of three ways, namely 
in fiber direction, in transverse direction and in shear. Therefore, we can suppose that at every 
current state of loading, the relative numbers of broken mesovolumes in cells for each type of 
failure are proportional to probabilities of failure rf. In the limiting case of the infinite number 
of mesovolumes in cell the proportionality becomes equality. So, we can interpret the numbers 
r = [rf, r 2 , r* 2 ] as relative counts of broken meso-elements of three types in each cell. 

Accumulation of damages in cells causes stiffness reduction of laminate and stress redis- 
tribution around the microcracks and between the cells and plies. This process is very compli- 
cated and modeling damage accumulation using micro-level consideration is, in general, not 
amenable to analytical studies. Probabilistic imitation approach [5] has lead to certain success 
in the area but it is computationally extensive and needs many simplifying assumptions to be 
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made as well. 


It is assumed in this paper that meso-element failure in certain direction causes deteriora- 
tion of properties in the same direction. Failure of r* fraction of all meso-elements in i-th direc- 
tion of k-th cell will result in reduction of cell elastic modulus in this direction. It is proposed 
that mathematical mean values of cell elastic moduli decrease with damage accumulation ac- 
cording to following hypothesis. 


E\ (r) = E? lo ( 1 - r k n (/) ) E\ (r) = E\ (1-4(0) 

Gn (0 = 4 0 (l-^ 2 (/)) ' vj 2 (r) = 4 0 ( 1 - 4 (0 ) 

but relative standard deviations remain constant 


eUo 




J D G*,(») _ J °12 0 

G 12 (0 G,2 o 


| 2 (') _ J V S 
4 0 


( 13 ) 


According to this algorithm damage accumulation in cell during laminate loading evokes 
the shift of cell moduli distributions towards zero direction and narrows these distributions pro- 
portionally to stiffness decrease. Changes in statistical distributions of stiffness resulting from 
damage accumulation according to the proposed algorithm of stiffness reduction are illustrated 
in Figure 2. 



figure 2. Stiffness distributions at various damage levels. 


This simple approach does not take into consideration micromechanical phenomena like 
stress redistribution around the single cracks or cracks interaction, but does take into account 
gradual stiffness reduction due to damage accumulation in cells and stress redistribution be- 
tween them. 


DISTRIBUTIONS OF CELL MATERIAL CHARACTERISTICS 

The basic feature of the proposed method is the consideration of stochastic nature of stress - 
strain field in the laminate due to stochastic stiffness and strength properties of laminate con- 
stituents. In addition, the concept of mesovolume as a building element for textile composites 
is introduced. Initial random properties of cells should be known to perform analysis. They can 
be obtained experimentally on small unidirectional specimens or can be calculated using some 
appropriate micromechanical models. In textile laminate, the scatter of angle of cell orientation 
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in ply should betaken into account. Experimental inspection of specimens of fabric laminates 
shows that scatter of orientation depends essentially on the weave type and manufacturing tech- 
nique. 

Experimental statistical data for Kevlar/epoxy unidirectional composite presented in ref- 
erence [16] are used in the following analysis. The random elastic and strength (critical strain 
bounds) characteristics are listed in Tables 1 and 2 (first line corresponds to zero deviation of 
angle). The standard deviations of random characteristics were calculated from their 95% con- 
fidential intervals given in [16], using well known statistical formulas. 

Effect of scatter of orientation on distribution of cell properties was analyzed using the 
Monte-Carlo method. First, properties of unidirectional material and angle orientation were 
chosen at random according to their distributions. Then cell properties were calculated 
using rota tional dependencies of the elastic and strength characteristics. Sets of several thousand 
samples were computed for each angle deviation value. Resulting mean values and standard 
deviations of variables are shown in Tables 1,2. Some distributions of material characteristics 
obtained at 5 deg orientation deviation are shown in Figures 3,4. 






figure 3. Calculated distributions of cell elastic characteristics (points) and their normal ap- 
proximations (lines) 

Analysis shows that scatter of orientation influences considerably the distributions of ma- 
terial properties. The strongest effect is observed for modulus £ t and Poisson’s ratio v 12 of the 

material. The increase of scatter of cell orientation causes decrease in tensile stiffness and slight 
increase in shear stiffness and failure strain of cell. In all cases, essential increase in deviation of 
material parameters is observed. Shape analysis of the resulting distributions shows that at 
larger deviations, some of them become significantly non-Gaussian. The reason for that is an 
extremal dependence of material properties on orientation angle around zero orientation. Fig- 
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table i. Mean Value and Deviation of Cell Elastic Characteristics for Different Orientation 
Scatter 


Angle St. 
Deviation 
deg 

GPa 

e 2 

GPa 

g 12 

GPa 

V 12 

0 

74.9/2.22 

4.65/0.16 

1.88/0.03 

0.35/0.031 

1 

74.0/2.48 

4.65/0.16 

1.88/0.03 

0.35/0.032 

2 

71.9/4.28 

4.65/0.16 

1.88/0.03 

0.36/0.039 

5 

62.0/12.9 

4.63/0.16 

1.91/0.05 

0.46/0.14 

10 

47.1/20.5 

4.60/0.16 

2.03/0.22 

0.71/0.37 


* mean / standard deviation 


table 2. Mean Value and Deviation of Cell Failure Strains for Different Orientation Scatter 


Angle St. 
Deviation 
deg 

eii 

% 

e" 

b ii 
% 

E 22 

% 

p" 

c 22 

% 

7'l2 ’ 7 "12 
% 

0 

1.71/.2 

-.47 8/. 024 

.283/.016 

- 1 .41/. 1 1 

2.561.19 

1 

1 . 11/2 

-. 478/.024 

.284/.016 

-1.41/.1 1 

2.55/.S2 

2 

1 . 12/2 

-.478/.024 

.283/.015 

-1 .41/. 1 1 

2.60/.76 

5 

1 . 12/2 

-. 482/.026 

.285/.016 

- 1 .42/. 1 1 

2.49/.78 

10 

1 . 16 / 2 ] 

-.494/.035 

.292/.021 

-1.45/. 13 

2.35/.92 


* mean / standard deviation 





ure 5 shows the variation of elastic characteristcs of unidirectional composite with orientation. 
Figure 6 represents the variation of tensile and compressive deformations at failure with angle, 
and Figure 7 - the variation of shear failure strain. These extremal dependencies cause the 
shape change of distributions. In case of low initial scatter the distributions may become one- 
sided. 




figure 4. Calculated distributions of cell failure strains (points) and their normal approxima- 
tions (lines). 




figure 5. Variation of cell elastic characteristics with fiber orientation. 




figure 6. Variation of cell failure strains in tension e' and compression e" with fiber orienta- 
tion. 

The performed calculation of third and fourth moments of distributions (i.e., coefficients 
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of skewness and kwitosis) and their comparison with literature data [17] have shown that for 
the material under consideration the 5 deg deviation of angle makes some of material distribu- 
tions to deviate significantly from normal distribution. Further development of method is re- 
quired to deal with these situations. 




figure 7 . Variation of cell failure strain in shear with orientation. 

ILLUSTRATIVE EXAMPLE 

Computer code is developed for calculating laminate response under the in-plane loading. 
Incremental loading is set in the program according to given derivatives 6, and equation (1) is 

integrated numerically. On each step of loading the damage functions in cells rf accumulated 
up to the previous step are used to calculate the current elastic properties of cells and laminate 
strain increments. It is assumed that complete failure of composite occurs when any of the cur- 
rent effective elastic moduli of laminate £,, E 2 , G n becomes equal to zero (in program - equal 
to small number). The program is written in “Mathematica” and handles general case of arbi- 
trary ply number, weave and arbitrary cell orientations in plies. Detailed information on dam- 
age functions and data on laminate strains and current effective properties are recorded to the 
data file on each loading step and plotted after completing the computation. 

The capability of program is illustrated on an example of Kevlar/epoxy laminate consisting 
of identical textile reinforced plies. Orientation of cells in the laminate is [±45] . Characteristics 
of Kevlar/epoxy composite with 2 deg orientation deviation (Tables 1,2) are used as cell mate- 
rial properties. The following values of weave parameter v were used in the analysis - 0.5, 0.33, 
0.25, and 0. 1 1, which correspond to plain weave, twill weave, 4 harness satin and 8 harness sat- 
in weave, respectively. 

Calculated stress-strain diagrams due to shear loading and cumulative damage functions 
of the composite are shown in Figure 8. Continuous lines, dashed lines and dash-dot lines in 
the following figures correspond to 1 1, 22 and 12 variables, respectively. Arrow points toward 
the lower values of v. Figure 9 shows the evolution of damages of different type in the lami- 
nate, obtained by averaging the damage functions over all cells. One can see that the more in- 
homogeneous distribution of cells between sublaminas is (lower v), the greater delay is 
observed for severe damage formation onset. The axial compressive deformations of laminate 
seen on the stress-strain diagrams (Figure 8) are due to damage induced anisotropy effect (see 
discussion in [15]). 





figure 8. Stress-strain diagrams and cumulative damage functions of textile laminate at shear. 



J 20 40 60 

x 12 , MPa 


figure 9. Average accumulation of damages of three types in textile laminate. 

CONCLUSIONS 

Damage evolution in te xtile laminates under the in-plane loading is discussed in this paper. 
A laminate model based on the special lamina subdivision is presented which predicts variation 
of properties with weave parameter. The damage accumulation in such a system is calculated 
based on the random function theory. Numerical algorithm and computer code for damage ac- 
cumulation and deformation history prediction are developed. The effect of orientation scatter 
on the statistical distributions of cell material properties is discussed. 
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